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ABSTRACT 

Near-surface coals in the Alberta plains occur in Upper 
Cretaceous and Tertiary sediments, the former represented by 
the Belly River Group, Horseshoe Canyon Formation and Wapiti 
Formation, and the latter by the Scollard Member of the 
Paskapoo Formation. These coal-bearing units contain 
Significant coal resources of varying quality and rank. 
Statistical techniques were used to document both the 
distribution of, and the interrelationships between, the 
components of proximate and ultimate analysis and calorific 
value determinations. 

The coals discussed range in rank from subbituminous C 
to high volatile bituminous C. These coals have variable ash 
content and low sulfur content. The relationships between 
calorific value (dry basis, kJ/kg) and ash (dry basis, %), 
and calorific value (moist mineral matter free, kJ/kg) and 
equilibrium moisture (MEQ, %) were determined to be 

cv(p) = 29262 - 286(Ash[D]) (+1642 kJ/kg) 
and CV(MMMF) = 31816 - 442(MEQ) (+1408 kJ/kg), 
respectively. 

A model for coal rank distribution in the Alberta 
plains, proposed in this study, took into consideration the 
maximum (paleo-) depth of coal seam burial. This maximum 
depth (DOB, m) was reconstructed on the basis of previously 
published graphs, relating the equilibrium moisture loss to 
depth of coal seam burial, uSing the relationship 


log[10]MEQ = 1.865 - 0.000416(DOB). 
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A least squares regression analysis of calorific value 
(MMMF, kJ/kg) on the calculated depth of burial (DOB, m) 
yielded the equation 

CV(MMMF) = 14748 + 6.25(DOB), 
indicating a coalification gradient, determined on the 
reconstructed overburden, of 6.25 kJ/kg per metre (0.82 
Btu/lb/ft) of depth (or overburden). 

Near-Surface coals in the Alberta plains increase in 
rank toward the west southwest (toward the foothills / 
mountains region) because of progressively greater amounts 
of overburden that existed in that direction during Tertiary 
time. Erosion has since removed between 900 and 1900 m (2953 
and 6234 ft) of sediment, with the greatest amount of 


removal occurring in a west southwest direction. 
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I. INTRODUCTION 


A. BACKGROUND 

In the course of its on-going investigation of the coal 
resources in the Alberta plains, the Geological Survey 
Department of the Alberta Research Council has accumulated a 
large number of proximate and ultimate analyses of coal 
samples. Between 1974 and 1980, 685 coal samples were 
analyzed and the results of these analyses were made 
available to the author for a research project (Appendix 1). 
The purpose of the project was two-fold. The first was to 
describe and summarize the distribution of each of the 
components of a proximate and ultimate analysis for each 
coal-bearing unit, using the recently acquired data. The 
last comprehensive attempt at this was by Campbell (1964), 
uSing data accumulated by the Alberta Research Council 
between 1944 and 1963. In addition to describing the 
distributions of individual components, component 
interrelationships were also to be investigated. The second 
purpose of the research project was to use the bed moisture 
content of the coal, which had previously been identified as 
possibly controlled by original depth of coal seam burial 
(Hacquebard, 1977), to reconstruct the maximum 
paleotopographic surface in the central plains area of 
Alberta. This surface can then be looked at as a possible 
controlling factor in determining the present day 


distribution of coal ranks. Coal geologists working in the 
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plains area of Alberta have generally not excepted this 
hypothesis, with these geologists tending to relate the coal 
rank variations observed in the plains area to forces 
accompanying the uplift of the Cordillera (Campbell, 1983; 


Campbell, 1972; Stansfield and Lang, 1944) 


B. REGIONAL GEOLOGY AND SAMPLE DISTRIBUTION 

The coal-bearing portion of Alberta can be divided into 
two geographic regions, the folded and faulted foothills and 
mountain region, and the flat-lying to gently dipping plains 
region (Fig. 1). Coals of significant thicknesses are known 
to occur in both areas in Lower Cretaceous, Upper Cretaceous 
and Tertiary rocks (Williams and Murphy, 1981; Yurko, 1975; 
Stansfield and Lang, 1944). 

Four prolific coal-bearing horizons are present within 
Upper Cretaceous and Tertiary sediments in the plains of 
Alberta: Foremost and Oldman Formations of the Belly River 
Group, Horseshoe Canyon Formation, Wapiti Formation, and 
Scollard Member of the Paskapoo Formation (Fig. 2). The 
outcrop extent of each of these units is shown in Figure 1. 
Except in areas of present-day active erosion (for example, 
the banks of major rivers), the entire plains area is 
covered by a thin layer of glacially-deposited Quaternary 
sediment. The coal-bearing strata dip regionally to the west 
southwest (ERCB, 1982), from near horizontal in the eastern 
portion of the plains area to over 8.6m per km in areas 


near the foothills (Williams and Burk, 1964). 
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Figure 1: Outcrop distribut; 

; ribution of Upper Cretaceo 
Tertiary coal-bearing units (from the Geolo a ae 
of Alberta). er orree 
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Figure 2: Upper Cretaceous and Tertiary Stratigraphic 
nomenclature of the central plains area. The 
Stratigraphic units noted were compiled following 
William and Burk (1964), Irish (1970) and Nurkowski 
and Rahmani (in press). 
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The Belly River Group', subdivided locally into the 
Foremost and Oldman Formations, and the Wapiti Formation are 
used as defined in Williams and Burk (1964), and the 
Scollard Member of the Paskapoo Formation is used as defined 
in Irish (1970). In this study the Horseshoe Canyon 
Formation encompasses the Horseshoe Canyon Formation as well 
as the Whitemud Formation of Gibson (1977) and Irish (1970). 
Use of this new definition of the Horseshoe Canyon Formation 
can be found in Nurkowski (1980) and Nurkowski and Rahmani 
(in press). 

The Belly River Group (Fig. 2), lowermost of the 
outcropping coal-bearing units, is sandstone dominated, 
containing coals, siltstones and shales. Individual coal 
seams attain thicknesses up to 3 m (10 ft), with cumulative 
coal thicknesses up to 6 m (20 ft; Holter et a]., 1978). 
Deposition of the sediments occurred in continental, marine 
and marginal marine environments (Williams and Burk, 1964). 

The Bearpaw Formation (Fig. 2) is shale-dominated and 
contains minor amounts of sandstone. Thin coal seams of 
considerable lateral continuity have been found in the upper 
part of this formation, but whether or not they belong 
stratigraphically and genetically to the overlying Horseshoe 
Canyon Formation is questionable. The Bearpaw Formation was 
deposited under marine and marginal marine conditions 


(Williams and Burk, 1964). 


‘also known as the Judith River Formation 
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The Horseshoe Canyon Formation (Fig. 2) is sandstone 
dominated, with coals, siltstones and shales. Coal seams up 
to 5m (16 ft) thick and areas with cumulative coal 
thicknesses up to 12 m (39 ft) occur within this formation 
(Holter et a]., 1976). Sediment deposition is considered to 
have occurred in both continental and marginal marine 
environments (Rahmani, 1981; Nurkowski and Rahmani, in 
press.). 

The Wapiti Formation (Fig. 2) occurs in the 
north-central’ portion of the province (Pig. 1)°and is 
defined aS containing sediments of the Belly River Group and 
the Horseshoe Canyon Formation in areas where the Bearpaw 
Formation is absent. The Wapiti Formation is sandstone 
dominated, containing coals, siltstones and shales. Coal 
seams up to 3 m (10 ft) thick have been encountered within 
this’ formation (Chu',-4978) . 

The Battle Formation (Fig. 2) is a wide spread unit 
deposited in a swampy to lacustrine environment (Binda, 
1970) and is devoid of coal. This formation is very useful 
for separating the lithologically similar upper Horseshoe 
Canyon Formation from the lower Paskapoo Formation. It is 
readily recognizable in outcrop by its medium to dark 
purplish-grey to lime-grey color (Gibson, 1977) and, in 
subsurface geophysical logs, by its very low resistivity and 
higher gamma ray readings. Due to its thinness (commonly 7 
to 11 m [22 to 36 ft]; Nurkowski and Rahmani, in press), the 


formation was not plotted in the outcrop map (Fig. 1). 
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The Scollard Member of the Paskapoo Formation (Fig. 2) 
is the uppermost of the major coal-bearing units and 
contains both the thickest and the most continuous coal 
seams. It is defined as the generally finer-grained section 
of sediments overlying the Battle Formation and containing 
the thick coals traditionally named the "Ardley" coals. Coal 
occurrences in Paskapoo Formation sediments other than in 
the Scollard Member are rare. Coal seams up to 5 m (16 ft) 
thick, as well as coal zones with cumulative coal 
thicknesses up to 24 m (79 ft) occur in this member 
(Nurkowski, in prep.). Paskapoo Formation sediments and 
coals are considered to have been deposited within a 
continental environment (Carrigy, 1971). 

Between 1974 and 1980 the Alberta Research Council 
Grilled 427 testholes into the above-mentioned coal-bearing 
units in order to evaluate plains area coal resources. These 
testholes are concentrated in areas where the coal-bearing 
units are shallow, near their outcropping edge (Fig. 1). In 
total, these drilling programs accumulated information on 
89.6 km (293,963 ft) of drilled sediment and 2.1 km 
(6,890 £t) of cored sediment. Coal samples were taken 
periodically when coal seams were encountered during 
conventional open hole drilling. These samples were washed, 
bagged and sent for analysis. Portions of coal seams 
collected from cores were also bagged and sent for analysis. 
Contamination of coal samples (especially cuttings) by 


up-hole coal seams was not considered a major problem. 
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Inspection of caliper seers on geophysical logs indicate 
that shallow coals (less than 500 m) tend not to cave or 
wash out, and that this property of coal seams would 
minimize contamination. 

In total, 685 proximate analyses, 621 calorific value 
determinations and 104 ultimate analyses were performed 


(Appendix 1). 


C. PREVIOUS WORK 

The first detailed account describing the composition 
and quality distribution of coals in Alberta was by 
Stansfield and Lang (1944), who used mine samples for 
analyses. Their study found that coal quality (increase in 
calorific value, decrease in moisture content) increased 
toward the mountain front. This increase was attributed to 
"mountain building pressures" and, to a lesser extent, to 
geologic age and present depth of the coal seam (Stansfield 
and Lang, 1944, p. 9). 


Campbell (1964) summarized coal composition and quality 


Gata for coal mines that had been, or were, operating in the 


plains area of Alberta. Campbell (1972) used coal samples 


from oil and gas wells to determine a vertical coalification 


gradient, which was found to be 5.34 kJ/kg (moist ash-free) 
per metre of depth (0.7 BTU/1b/ft). In explaining the 


increase in coal rank in a direction toward the mountain 


front, Campbell (1983) supported the views of Stansfield and 
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Cordillera ... was accomplished by a continous horizontal 
force ... that could metamorphose the coal just as 
effectively as overburden weight." 

Steiner et al]. (1972), in part using data from Campbell 
(1964), also discussed composition and quality variation of 
coals from the plains area of Alberta. As in Stansfield and 
Lang (1944), moisture content was found to decrease, and 
calorific value to increase, towards the foothills and 
-mountain region. This study also noted that ash content 
increased to the southwest, an effect that was attributed to 
contamination from wind-blown volcanic ash from the southern 
Cordillera. 

Previous work on coal metamorphism, or coalification, 
has been well summarized by, among others, Bustin et al. 
(1983), Berkowitz (1979) and Teichmuller and Teichmuller 
(1968), and this study relied heavily on these references. 
Coalification results from the interaction of three 
controlling parameters: temperature, pressure and time. 
IncreaSing temperature is considered to be responsible for 
chemical changes associated with increasing rank (Berkowitz, 
1979: Teichmuller and Teichmuller, 1968) and increasing 
pressure is associated with physical-structural changes 
(Teichmuller and Teichmuller, 1968) that accompany 
increasing rank. The third parameter, time, Speeeee 
coalification in that a higher degree of coalification can 
be achieved during a shorter period of time if the coal was 


subjected to a higher than usual temperature. Relationships 
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between time and temperature have been developed by Karweil 
(1956) and are discussed in more detail in Hacquebard (1977) 
and Berkowitz (1979). These relationships were applied 
exenSively in Hacquebard's (1977) study of the coalification 
history of Mannville coals. Here, Hacquebard used bed 
moisture contents of near-surface coals from Steiner ef al. 
(1972) and the vitrinite reflectance of Mannville Formation 
coals (Lower Cretaceous) to estimate maximum burial depths 
and paleotemperatures for Mannville Formation sediments. 
Because vitrinite reflectances, which are indicators of 
thermal maturation (Teichmuller and Teichmuller, 1968), were 
not available, application of Karweil's (1956) model to 
determine the thermal history of the coals discussed in this 
study was ret possible. In order to obtain a rough estimate 
of the temperature required to mature the coals discussed in 
this study, the maximum vitrinite reflectance expected was 
estimated from a relationship of vitrinite reflectance 
versus carbon content from Berkowitz (1979). This value of 
vitrinite reflectance could then be compared to a table of 
minimum paleotemperatures required to mature coals to a 
Particular reflectance value (Bostick, 1979). 

Coal seam nomenclature, depositional settings, and 
Gescription of overlying and underlying lithology for 
coal-bearing horizons can be found in Rahmani (1981) for the 
lower Horseshoe Canyon Formation, Nurkowski (1980) and 
Nurkowski and Rahmani (in press) for the upper Horseshoe 


Canyon Formation, Chu (1978) for the Wapiti Formation, 
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Holter et al. (1978) for the Belly River Group, Gibson 
(1977) for the Horseshoe Canyon Formation and the Scollard 


Member, and Holter et al]. (1975) for the Scollard Member. 


D. COMPONENTS OF A COAL ANALYSIS 

A complete coal analysis generally consists of a 
proximate analysis, a calorific value determination and an 
ultimate analysis. 

Proximate analyses involves the determination of 
moisture, ash, volatile matter and fixed carbon, with each 
component expressed aS a percentage of the sample. These 
components can be expressed or reported in a number of ways, 
depending on the moisture value used. Samples expressed to 
an as received basis use the moisture value of the coal as 
it was submitted to the laboratory. The as determined basis 
of reporting uses the moisture value obtained during the 
analysis. The difference between these moisture values 
occurs because coal samples are generally dried before they 
are analyzed. Both as received and as determined moisture 
contents are provided in Appendix 1. Dry basis can be 
calculated from either the as determined or the as received 
basis by removing the moisture component and calculating the 
remaining components to total 100 percent. Dry ash free 
basis of reporting removes both the moisture aia Sen of the 
coal, and the remaining components, volatile matter and 
fixed carbon, are recalculated to total 100 percent. The 


method used to calculate the dry basis and dry ash free 
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basis is further explained in Appendix 2. 

Calor.fic values are reported in Appendix 1 in the 
metric form, kJ/kg. To translate these units to BTU's/1b, 
multiply the kJ/kg value by a conversion factor of 0.43. For 
comparison, calorific values can be calculated to both a dry 
and a ary ash free basis following methods outlined in 
Appendix 2. 

Ultimate analysis of coal includes determination of 
carbon, nitrogen, hydrogen and sulfur. Oxygen, also part of 
an ultimate analysis, is determined by difference. When 
these analysis are available, they are reported on a dry 
basis in Appendix 1. The oxygen and hydrogen contents 
reported in an ultimate analysis do not include the oxygen 
and hydrogen of the sample moisture. 

Equilibrium moisture or bed moisture values 
(Appendix 1) represent the moisture content of the coal in 
Situ (i.e. as it occurs in the seam). All components of a 
proximate and ultimate analysis, and calorific value, can be 
recalculated to a moist value (equilibrium moisture basis) 
following the method described in Appendix 2. 

In order to use the ASTM' rank classification system of 
coal (Fig. 3), it is important to distinguish between ash in 
coal, as determined by the test for ash, and mineral matter 
in coal. The ASTM classification system Saeiliiag Fee 
calorific value, volatile matter and fixed carbon be 
expressed on a mineral matter free’ 


‘American Society for Testing and Materials. 
2Calorific value must be expressed on a moist mineral matter 
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Calorific Value (kj/kg) 


Rank Equal to or Less Than 
Greater Than 


High Vol. Bituminous A 32564 * = 

High Vol. Bituminous B 30238 32564 

High Vol. Bituminous C 26749 30238 
24423 26749** 


Subbituminous A 24423 26749 
Subbituminous B 22097 24423 
Subbituminous C 19306 22097 
Lignite A 14654 19306 
Lignite B - 14654 


*Also Based on Fixed Carbon and Volatile Matter 
**If Agglomerating 


Figure 3: ASTM classification of coals by rank. 
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basis. During conventional ashing of coal, dehydration (e.g. 
FeSO,.nH20 to FeSO,), decomposition (e.g. CaCO; to CaO), and 
oxidation (e.g. FeS to FeSO,) occur, as well as partial loss 
of volatile constituents, particularly Hg, K, Na, Cl, P and 
S (Berkowitz, 1979, p. 49). The percent ash value obtained 
for the sample thus does not reflect the true mineral matter 
content in the sample. Stansfield and Lang (1944) claim 
that, on the average, ten parts of mineral matter in coal 
leave only nine parts ash. The Parr Formulae (Appendix 3) 
allow for translation of the analytical values to a mineral 
matter free basis. 

In Appendix 1 and in all descriptive and statistical 
sections that follow, dry ash free basis and mineral matter 
free basis of reporting have been calculated when the 
percent ash value (dry basis) was less than 25 percent. As 
the ash content increases, the volatile matter content, when 
compared to the fixed carbon content, also increases, and to 
avoid this an ash cutoff was required. This will be 
discussed further in Chapter II. 

The coal samples collected in 1974 and between 1976 and 
1980 were analysed by Chemical and Geological Laboratories, 
Edmonton, and those taken in 1975, by Cyclone Engineering, 
Edmonton. A number of samples from the above mentioned years 
were also analysed by the Alberta Research Council. No 
attempt was made to check the repeatability or accuracy of 
the analytical results from the various labs. Personal 


7(cont'd)free basis and volatile matter and fixed carbon 
must be expressed on a dry mineral matter free basis. 
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communication (1981) with each of the labs indicated that 


ASTM procedures were followed in all cases. 


E. METHOD OF STUDY 

In order to determine the coal-bearing formation from 
which each of the coal samples was collected, previously 
published Alberta Research Council reports were searched and 
coal seam correlations made in these reports were noted. In 
general, testholes drilled during 1974 evaluated the 
Scollard Member, during 1975 the Horseshoe Canyon Formation, 
during 1976 the Belly River Group, during 1977 the Wapiti 
Formation, and during 1978, 1979 and 1980 the Horseshoe 
Canyon Formation. 

In 1980 and 1981 the results of the 685 coal analyses 
(Appendix 1) were placed in the Alberta Research Council's 
Digital Vax 11/780 computer. A computer data base management 
system known as DATATRIEVE was used to retrieve data 
selectively and recalculate them to consistent bases. In 
order to compute both descriptive statistics on analytically 
determined components and component interrelationships, 
output from DATATRIEVE was transferred to the University of 
Nibebta’s Amdahl 470V/8 computer. Here, the MIDAS statistics 
package was used to generate the information contained in 


the tables, histograms and cross-plots that follow. 
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Pi PRESENTATION AND DISCUSSION OF DATA 
The analytical techniques used to determine moisture, 
ash, volatile matter and calorific value will be briefly 
described here. For a more detailed discussion, the reader 


is referred to Berkowitz (1979) and other coal handbooks. 


A. DEPTH 

The coals analyzed in this study were obtained from 
Grilling programs. The depths, summarized for each 
coal-bearing unit in Table 1, represent the top, or 
Shallowest depth of the sampling interval. As will be 
discussed later, quality characteristics of coals are 


related to both depth and geographic location of the sample. 


Table 1 - Summary of sample depth (in metres). 


Coal-bearing Unit n min max mean 
Scollard Member 227 20 275 168 
Horseshoe Canyon Fm 365 18 319 114 
Belly River Group 41 54 250 146 
Wapiti Formation 49 15 217 102 
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B. PROXIMATE ANALYSIS 


Moisture 

The most important moisture value in a coal sample is 
the equilibrium moisture, as it is taken to represent the jn 
situ moisture of the coal. Berkowitz (1979, p. 31) describes 
equilibrium moisture as 

",..the amount of water that a coal will hold when 

fully saturated at nearly 100% relative humidity (as 

in an undisturbed coal seam) and thus reflects the 

total pore volume of the coal accessible to 

moisture.” 
Equilibrium moisture content is determined following the 
procedure outlined in ASTM 1492. This requires that a sample 
of coal be brought to 97 to 98% relative humidity (at 30 
+0.2 degrees C), weighed, dried, and weighed again to 
determine the amount of moisture loss. 

This study, as well as previous research (Dyck et al., 
1980), noted an inverse relationship between equilibrium 
moisture content and ash content (discussed later). For this 
reason the descriptive statistics of equilibrium moisture 
content (Table 2) only consider coal samples containing less 
than 25% ash (D)'. The histograms of equilibrium moisture 
distribution, however, show all data (Fig. 4). | 

It should be noted that in some cases the mean of 


equilibrium moisture does not occupy the modal class of the 
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Table 2 - Summary of equilibrium moisture (in percent) for 
coals with less than 25% ash (D). 


Coal-bearing Unit n min max mean Sid. 
Scollard Member 72 See) iLo'e6 14.4 Ie6 
Horseshoe Canyon Fm 160 1253 24.7 1956 Zee: 
Belly River Group 16 8.5 26S 14.0 6:3 
Wapiti Formation 0 n/a n/a n/a n/a 


respective histogram (Fig. 4). This is because equilibrium 
moisture values are not normally distributed in the 


Statistical sense. 


Ash 

Ash is the inorganic residue remaining after coal 
combustion (Stansfield and Lang, 1944), and is determined 
following the procedures outlined in ASTM 3174. A weighed 
Sample of dried coal is combusted in excess air at a 
temperature of 725 (+25) degrees C, and the residue is 
weighed. As previously discussed, the residue remaining does 
not reflect the total mineral matter content of the coal. 
The latter is calculated using the Parr Formulae 
(Appendix 3). 

As seen from the histograms of ash distribution 
(Fig. 5), the ash (D) content of coal samples varies 
considerably. Samples containing more than 25 percent ash 
are probably argillaceous coals, partings or carbonaceous 


shales. After all such samples were removed, descriptive 
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statistics were calculated for all coal-bearing units 


(Table 13) 9 


Table 3 - Summary of ash (dry basis, in percent) for coals 


wi tielessiit hame25% ashe( pir, 


Coal-bearing Unit n min ma x mean Sid. 
Scollard Member 128 Shere 24.9 Seis 4.9 
Horseshoe Canyon Fm 205 he 24.9 IGS, 6. | 
Belly River Group 32 oe 24.7 1666 org 
Wapiti Formation 47 4.5 24.6 1185 4.5 


Volatile Matter 

Volatule matter “is the portion of the coaluthat is 
driven off during heaentla) and includes such substances as 
CO, COZ) CHS. tars and.other hydtogarbons. It fusadetermined 
following procedure ASTM 3175. A weighed sample of coal is 
heated in a covered crucible to 950 (+20) degrees C for 7 
minutes, with care taken to prevent oxidation of the sample. 
The remaining sample is weighed and the weight loss, minus 
the weight of sample moisture, determines the volatile 


matter constituent of the coal. 
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Due to the wide distribution of ash values (Fig. 5), 
volatile matter and fixed carbon values should be expressed 
on a ary ash free basis for comparison. However, when 
calculating high ash samples to this basis, the ratio of 
these components to each other can become distorted, in that 
in progressively higher ash samples, the volatile matter 
content increases and the fixed carbon content decreases 


CRabler4)\ 


Table 4 - Variation of the mean value of volatile matter 
content with respect to increasing ash content (all values 


in percent). 


Ash Volatile Matter Ash Volatile Matter 


(dry) (dry ash free) (dry) (dry ash free) 


Ova ne 42.4 50 = 60 49.3 
Dete2'5 43. | GOR 70 3 One 
1.03 ES 2'0 42.7 Oe = 80 62.4 
202-93 0 45.1 BO. = 90 76.4 
30 - 40 47.3 907 <2 100 82.5 
a0. -" 50 47.1 


By comparing the volatile matter content at the 0 - 10%, 
10 - 20% and 0 - 25% ash levels (Table 4), it can be noted 


that below 25% ash, the volatile matter increase is 
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relatively low, compared to coals with ash contents higher 
than 25%. The increase in volatile matter content is due to 
chemically bound water in the mineral matter portion of the 
Sample being driven off during the test for volatile matter. 
Clay dehydration curves (Grim, 1953) indicate that most 
chemically bound water is driven off at temperatures between 
400 and 600 degrees C, a temperature considerably higher 
than encountered during the test for moisture (107 degrees 
C). Taking this into consideration, histograms and 
descriptive statistics for volatile matter (DAF)' are 
prepared for samples containing less than 25 percent ash 


(Fig. 6 iand-Table 5). 


Table 5 - Summary of volatile matter (dry ash free basis, in 


percent) for coals with less than 25% ash (D). 


Coal-bearing Unit n min max mean Sree 
Scollard Member 124 34,1 58.3 40.6 Stok 
Horseshoe Canyon Fm 137 34.6 65.4 44.0 31.0 
Belly River Group S2 S9e2 48.4 44.6 2n3 
Wapiti Formation 47 30.9 49.2 44.2 Pie: 
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Figure 6: Distribution of volatile matter (dry ash free 
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_ Fixed Carbon 
The term fixed carbon is used to describe that portion 
of coal remaining after moisture, volatile matter and ash 


have been removed. It is calculated by 


“FC = 100 - (%moisture + %ash + %VM). 


As discussed previously (and in Appendix 2), values for 
volatile matter and fixed carbon, expressed on a dry ash 
free basis or a ary mineral matter free basis, total 100 
percent. The descriptive statistics for fixed carbon (Table 
6), therefore, are related to those of volatile matter. Only 
coal samples containing less than 25% ash (D) are 


considered. 


Table 6 - Summary of fixed carbon (dry ash free basis, in 


percent) for coals with less than 25% ash (D). 


Coal-bearing Unit n min max mean s.d. 
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C. CALORIFIC VALUE 

Calorific value of a coal sample is determined 
following the procedure outlined in ASTM 2015. This requires 
that a weighed sample of coal be combusted in a calorimeter, 
with the calorific value of the coal sample being equated to 
the resulting rise in water temperature. Values are 
expressed in kilojoules per kilogram of coal. 

In order to use the ASTM classification of low rank 
coals (Fig. 3), calorific values must be expressed on a 
moist mineral matter free basis. Descriptive statistics for 


as well as for calorific values 


these calorific values, 


expressed on a ary ash free basis, are presented in Table 7. 


Table 7 - Summary of calorific value (in kJ/kg) 


with less 


thaniezoe ash .(D) 


for coals 


Coal-bearing Unit base n min ma x mean s.d. 
Scollard Member MMMF V2) BL 2OSS 22) O50. 125,479 796 
DAF “126 206@84. S317599 29,926 636 
Horseshoe Canyon Fm MMMF 155 20,204 27,130 23,048 1245 
DABS 9) ae OPES op eS O29 30 868 
Belly River Group MMMF lose 8207108 §—29--3565.26, 065. 2974 
DAF Sl 205 4609 32, O94 29, cco, 1 a9 
Wapiti Formation MMMF 0 n/a n/a n/a n/a 
DAF £99 26,549 32,594 720, 540 849 
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The distribution of calorific values, expressed to the above 
mentioned bases, are provided in Figures 7 and 8. In these 
histograms and descriptive statistics, only coal samples 
containing less than 25 percent ash were considered. 

As can be noted from Table 7, the coals discussed range 
in rank from subbituminous C to high volatile bituminous C. 
The mean carbon (DAF) content of these coals varies from 
73.0%, DOr Subbituminous,C coals to 7639% for high volatile 
bituminous C coals. In order to determine an estimate of the 
vitrinite reflectance values that would be expected for 
coals of these ranks, the values of carbon (DAF) content 
calculated were applied to a vitrinite reflectance versus 
carbon content relationship provided by Berkowitz (1979). 
This relationship indicates that the vitrinite reflectance 


Values for these coals should range fromM@Om4iyto 0.5% 
D. ULTIMATE ANALYSIS 


Susrur 

Because of the association of burning high sulfur coals 
and environmental pollution, there is considerable interest 
in the sulfur content of coals. The totald sulfur content, 
both organic and inorganic, was determined for a significant 
proportion of the samples. These contents are soneen on a 
dry basis in Appendix 1 and Table 8. 

The sulfur (D) distributions for the Scollard Member, 
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Table 8 - Summary of sulfur (dry basis, in 


coals. 
Coal-bearing Unit n min 
Scollard Member 199 ets 
Horseshoe Canyon Fm 358 O20 
Belly River Group 41 Oss 
Wapiti Formation 49 On" 


95% confidence interval on the 


Scollard Member 0°36 
Horseshoe Canyon Fm 0.44 
Belly River Group 0.60 
Wapiti Formation 0.29 


positively skewed with a positive kurtosis 


reason the Belly River Group does not have 
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distribution may be a function of the number of samples. 


More normal distributions occur if the sulfur percentages 


are transformed to their log[10] values. The minimums, 


maximums, means, and 95 percent confidence 


interval about 


the sample means, presented in Table 8, are calculated from 


the transformed data. 


Carbon, Hydrogen, Nitrogen and Oxygen 


The mean value of carbon, hydrogen, nitrogen and 


oxygen, on both a ary and dry ash free basis, are presented 


for each of the coal-bearing units in Table 9. The number of 
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analyses performed for these components was too few to 
permit the calculation of meaningful descriptive and 


comparative statistics. 


Table 9 - Mean value of ultimate analysis components (in 


Coal-bearing Unit base C H N 
Scollard Member DRY 6207 36 0) 
DAF BS 06 4.6 22 

n an Sy 26 

Horseshoe Canyon Fm DRY SO.55 4.1 2 
DAF 74.5 Om! AS 

n 53 53 39 

Belly River Group DRY 5 laut i 5) 
DAF 74.7 ree oo 

n 10 10 4 

Wapiti Formation DRY 63.8 4.3 aC 
DAF uel 38 4.8 a5 

n 4. & 4 
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E. COMPARISON OF RESULTS FROM CUTTINGS AND CORE SAMPLES 

The analyses used in this study were of samples 
collected from both cuttings and cores from conventionally 
drilled testholes. Cuttings were subjected to a 1.8 specific 
gravity separation with the less than 1.8 Specific gravity 
fraction analysed. This was done to remove the non-coal 
contamination portion of the sample. Core samples were 
analysed as sampled. A statistical bias limiting meaningful 
comparisons between the two sample types occurs when the 
above is taken into consideration. The ash content is 
artificially lowered in the cuttings samples (cuttings mean 
= 22.1%, core mean = 34.3 %); thus all parameters relating 
to ash content, particularly calorific value (as will be 
discussed later), will also be affected in the comparison. 
Rank determination, however, will be unmarred by this 
sampling method, as the calorific value is calculated to a 
mineral matter free basis. 

A statistical comparison of the two sampling methods 
will not be attempted here. The interested reader is 
referred to Dyck et a]. (1980) for their discussion of a 


Similar comparison. 
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III. DATA INTERRELATIONSHIPS 


A. CALORIFIC VALUE VERSUS ASH 

A close inverse relationship exists between the 
calorific value of a coal sample and its corresponding ash 
content. A crossplot of these variables, both expressed on a 
dry basis, is provided in Figure 10. All data with ash (D) 
content less than 25% are plotted, regardless of 
coal-bearing unit. A least Square regression analysis of 
calorific value (D, kJ/kg) on ash (D, %) yields the 


equation' 
CV(D)h=529269 = 286(ASH[ Dido # (le 


The correlation coefficient (R) between these variables is 
-0.89 (with R? = 0.78), and the 95 percent confidence 
interval on calorific value estimation, following the method 
described in, Ezekiel and Pox (01959) F359). Ziad depend S 
+: 1642 kJ/kg,(706 Btu/1b). 

The calorific value data shown in the crossplot 
(Fig. 10), and used in calculating the relationship between 
calorific value and ash, are from the entire plains area; 
thus the constant shown in Equation (1), representing the 
calorific value at a theoretically occurring zero ash value, 
would represent an average calorific value for near surface 
coals throughout the area sampled. 


'For Equation 1, the 95% confidence interval on the slope is 
+14.6, and on the intercept is +227.0. 


34 


qaate 2ntels stitaa sts, ter? sce. ,cek bis set 7 
eit gniddseoygs: ,{1)'\ wocteepe int mode (Im8FEROD 5 
euler ige ote: Pniziueos yi lasti¢siosds: ne? vig a8 tee rks 


eoiixue s65n toi aglav Si tivo ofsa apET OVA Ae. pice am to 


sefielbiev aacd2) neawitea (2), taeiatiteas weal aNOD a 
sonebitnoo inseise@e say Baw ,T8Y 0 = 'R gow) ca: 
my 


oi poiwollo® ,nanreniaes sulsy a tijokes Xe sem 
r, 


, . I j 
sa Otdanoljele: sid eniseiuates pe peat bne ADR 


As “e 1s fo Wier, a a 
. = 220, en eve oe 
7 — e ¢ 
~ i’ 6 Ley Foxe 7 
| = 7 fen = 
: > 
| SUERAMGI TALIA, pao Srey oe 
” ie 
x : " ¥ _ [ ro a 
a f - bey, see eee. 
Ben BUSHEY TutAy OTERO 
Bei ’ ~~ > =e ah 4 : acs rt Es a. 
ifevted e@tarrxs ¢ ia enoise ley .oetavags 29080) ey 
re | : oe 
[Bnogessaioo eat bas alqmag Dao & To Susy, Stee eres 
2 ot 
as " siatel BY te er yr) 0 4 ie ete) P| & 
if a ¢ : 
de C8 « a ak Jig ¥ ‘ '. 35 vor ¢ ~a 4 2 
pes oI By - CEsolg 6 #65 Aa evel 3 
i 
: e . ipat SaeRave, Faz /* 300 o 
iy 2) Aes. eo 2A iy guts 
* 
: Ge |H@A) ans, ~* Sekes x {GW i 
A / 
3 " ? 
aay oa 


' 
rae 


“ - Wee: 
aa8t} eat Bia! fos tena Abd apie: 
@ ok. : 1 

aes Nose 20°) ae 3088 


” 

* 

. 

4 
fa. 
o~-“- 

Pr 

| 

os 

a 


Jobqeéors, eae nd dead B25 wiley our a4 a oe: 
ms 1 


f ' ~~ - 


; an 
ae: che 
ie : age at at 
a 
iei¢eed 
™ i 
. ay 


eae 


* - ' a) fi hy 
“ i oS i » 
ae ve 
a 
4 ja 5 a 
" 


ide i ll Bar8) ri #491 * te | 


iia 


ir e vo 


: * bE 0. €$S:3 : ok 349 DI STN s 


b onpty: 0s ‘ances Diibigciaae (ae a a 
12 sgole- sit no tawssn ni sonebitago Ree AF .', 1 ORS 


° 
29,000 5 6 
° ° e. ° Single Data Points 
° 
° 
0 0 ° e Multiple Data Points 
28,00043° °° oe ° 
° 
° 96 Pole coco ° Oo N = 406 
B Gte. os C0 00 0 oO ° 
27,000 Gost e Perens ° 
° 
sek ° rae @0 @ on 0° ° 
~ ° ° 0 00 ° 
= ° °° (ooo) 
= : eo 00 ° ° 
ray 26,000 . pers oc} (oe oe) ° P 
ew eo Oo co @ ° ° 
@ (ey fete) co ©0@08@ 00 
=) ° % oe oe ° 
i 25,000 OnsmeOu HO ° © 0 o8®\_ce~ 8 ooo ° 
> ° r) ° ° ° ° 
° ° ® fore) Ie) QO @0o ° 
oO ° oO ° ;. ° 0° ° 
aS : ° 
— ° bo 3 0° i ° *: 
~ 24000 oe) CODON @ © 0 O ° 
Oo ¥ ° ° @o oo 
r @®0° 0 000 OY 0 CO © © 
e o 
O 4 a ae ee 
23,000 oG Meger Coote 
° (oe) © @O\.0 60 
oO oo 
° 
2) 8) : 
22,000 a mto te teas ao: 
° ° 
° ° 
° 
Orn @o 
21,000 ° 
° 
5 10 15 20 25 


Ash (D), % 


Figure 10: Crossplot of calorific value (dry basis) 
versus ash (dry basis). 
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The slope component of Equation (1) indicates that dry 
calorific value decreases by 286 kJ/kg (123 Btu/lb) per 
percent increase in ash (D) content. 

Table 10 shows the results of both correlation 
coefficient and least squares regression analysis 
calculations for each of the coal-bearing units. Of 
Particular interest is the lower correlation coefficient 
between ash and calorific value for the Belly River Group 
(R = -0.65) compared to the remaining coal-bearing units 
(with R = -0.87 to -0.94). The reason this occurs may be 
that the Belly River Group coals were sampled over a 
considerable area in an east - west direction, from township 
3, range 2 W4M to township 9, range 22 W4M (Fig. 2), whereas 
Horseshoe Canyon Formation sndvseottaca Member coals were 
Sampled in a narrow north - south band following the 
outcropping edge of the formation. The east - west variation 
in the Belly River Group samples introduces a change in 
calorific value because of increasing rank (hence an 
increasing degree of maturation) in a westerly direction. 
Therefore, locations of Belly River Group samples affect the 
ash / calorific value relationships more than for other 


coal-bearing units. 


B. EQUILIBRIUM MOISTURE VERSUS ASH 
The calculated correlation coefficient (R) between 
equilibrium moisture and ash (D) content is -0.38, which 


indicates that samples with higher ash have lower 
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Table 10 - Equations relating calorific value (dry basis, in 
kJ/kg) and ash (dry basis, in percent) for coal samples 
containing less than 25% ash. 


Coal-bearing Unit n R Equation 

Scollard Member (2650594 CV=29955-301(ASH) 
Horseshoe Canyon Fm 197 -0.94 CV=29668-318(ASH) 
Belly River Group 62.4 =00 65 CV=27720-203(ASH) 
Wapiti Formation ee Or Bi CV=28375-286(ASH) 


equilibrium moisture. This value of correlation coefficient 
(R = 0.38) is statistically significant with 419 samples 
used in its calculation. 

Dyck et a]. (1980) calculated a correlation coefficient 
of -0.90 for equilibrium moisture and ash (D) content, using 
1176 samples of Saskatchewan lignites. This is considerably 
higher than the correlation coefficient calculated in this 


study. A possible reason for this will be discussed later. 


C. CALORIFIC VALUE VERSUS EQUILIBRIUM MOISTURE 

The calorific value of a coal sample, expressed on a 
moist mineral matter free basis is affected by the 
equilibrium moisture of the sample, thus correlation between 
calorific value (MMMF)' and equilibrium moisture is to be 
expected. A correlation coefficient (R) of -0.91 (R? = 0.83) 
was calculated for these variables for all data regardless 


of coal-bearing unit. A cross-plot of these variables is 


‘moist mineral matter free basis. 
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presented in Figure 11. A least squares regression analysis 
of calorific value (MMMF, kJ/kg) on equilibrium moisture 


(MEQ, %) yields the equation’ 
CV(MMMF) = 31816 - 442(MEQ) (2). 


The 95 percent confidence interval on calorific value 
estimation is + 1408 kJ/kg (605 Btu/l1b). 
| This relationship (as shown in Fig. 11) indicates the 

Significant degree to which moisture content and coal rank 
are inter-dependent in low rank coals. Teichmuller et a]. 
(1968, p. 248) similarly state that 

".,..according to recent work the moisture content 

and the moisture-holding capacity of ash-free coal, 

as well as its calorific value which depends to a 

large degree upon the moisture content, remain the 

best rank indicators up to the stage of 

high-volatile B and A coals". 
The slope component of Equation (2) indicates that calorific 
value (MMMF) decreases by 442 kJ/kg (190 Btu/lb) per percent 


increase in moisture content. 


D. VOLATILE MATTER AND FIXED CARBON VERSUS CALORIFIC VALUE 
The correlation coefficient (R) between volatile matter 

(DAF) content and calorific value (DAF) is -0.31, indicating 

that an increase in calorific value is associated with a 


‘For Equation 2, the 95% confidence interval on the slope is 
+25.2, and on the intercept is +456.8. 
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Figure 11: Crossplot of calorific value (moist mineral 
matter free basis) versus equilibrium moisture. 
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decrease in volatile mar demieon: ene (or conversely with an 
increase in fixed carbon content). 

These relationships have been known for some time and 
are discussed in detail by Teichmuller et a]. (1968) and 
Berkowitz (1979). The values of correlation coefficient 


(R = +0:31) are statistically significant with n = 402. 
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IV. ANALYSIS OF RESULTS AND PREDICTIVE MODELS 

Hilt's rule indicates that in undisturbed strata, the 
rank of coal increases with increasing depth of burial, and 
this relationship has been discussed by many authors 
(Berkowitz, 1979; Hacquebard, 1977; Stach et al]., 1975; 
Teichmuller et a/]., 1968). In undisturbed strata an 
increasing depth of burial facilitates coalification in two 
ways, by increasSing temperature, at a rate equal to the 
geothermal gradient, and by compaction, through increasing 
overburden pressure. The increaSing temperature is 
considered to be responsible for chemical maturation 
(Berkowitz, 1979) and the degree of maturation can be 
determined by meaSuring the vitrinite reflectance. 
Reflectance data, unfortunately, were not available for the 
coals discussed in this study. Thus comparing the regional 
coal maturation levels attained, and determining the maximum 
paleotemperatures to which the coals were subjected (through 
the use of Karweil's (1956) relationships), were not 
possible. 

A very rough estimate of the minimum temperature 
required to mature the highest rank coals in this study can 
be obtained from Bostick (1979). Bostick (1979) provides the 
minimum paleotemperatures required to mature coals to 
various degrees of reflectance values (starting ae 0.48% 
reflectance). As was earlier determined, the coals in this 
study could have vitrinite reflectance values as high as 


0.51%. This value of reflectance corresponds to a minimum 
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paleotemperature of 107 degrees C. It should be noted that 
Bostick (1979) suggests the table should be used for coals 
Mesozoic and older in age. The high volatile bitminous C 
coals in this study vary in age from late Mesozoic to early 
Cainozoic (Campanian to Paleocene), so the minimum 
temperature estimated should be used with caution. The 
estimate of minimum paleotemperature from Bostick (1979), 
however, iS similar to the calculated paleotemperature using 


the equation 
Ths 302.97(logLiOlRe) etsie@z.233 


from Price (1983), which is 99 degrees C. Contrary to many 
coalification models, Price (1983) disregards time as an 
important function during coalification, and the difference 
between the two temperatures (107 versus 99 degrees C) may 
be related to this matter. 

The previous discussions on estimating 
paleotemperatures are based on estimates of vitrinite 
reflectance values from carbon (DAF) contents. Because 
actual reflectance values were not available, this study did 
not investigate the chemical changes that accompany 
increases in temperature. However, the data are available 
for studying the increases in rank due to cempaction, which, 
in a tectonically quiescent setting such as the plains area 
of Alberta, should be a function of the depth of coal seam 


burial. 
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The relationship between increasing rock density (or 
decreasing pore volume) and increasing depth of burial has 
been noted by many authors (Magara, 1976; Till, 1974; 
Teichmuller and Teichmuller, 1968) and similar compaction 
phenomena should occur in coals. Berkowitz (1979) associates 
equilibrium moisture content with porosity (pore volume) and 
Hacquebard (1977) relates decreases in equilibrium moisture 
content to increases in original depth of burial’ 

(Fig. 12)?. Using this relationship, and equilibrium 
moisture data from Steiner et a]. (1972), Hacquebard (1977) 
calculated the thickness of removed sediment over oil and 
gas fields in Alberta. This study, using a regionally 
distributed and recently collected data base, will try to 
oe the area over which the thickness of eroded sediment 
can be calculated, and relate this original thickness of 
sediment to the rank variation patterns noted in the plains 
area of Alberta today. 

Hacquebard's (1977) graph (Fig. 12) indicates a 
logarithmic decrease in bed moisture (equilibrium moisture 
[MEQ, %]) with respect to an increase in depth of burial 
(DOB, m). By eStimating values from both axes of the graph, 


this graph can be translated to the formula 


Tog l0i}MEO = 17865 — 0.000416(DOB) (3). 


‘Also see Teichmuller and Teichmuller for relationships of 
equilibrium moisture, depth of burial and rank. 

*The data used by Hacquebard (1977) to construct the depth 
of burial / moisture relationship (provided in this study as 
Figure 12) was derived from European research. 
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Depth of Burial in Feet 
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Figure 12: Relationship of moisture content in coal and 
its original depth of burial (from Hacquebard, 1977). 
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It should be noted that eae eos rank divisions indicated in 
Figure 12 are slightly different from those encountered in 
this study. The relationship between equilibrium moisture 
and depth of burial for coals in the Alberta plains coal 
will be discussed later. 

Using Equation 3, the original depth of burial was 
calculated for all coal samples with ash contents less than 
25% for which equilibrium notecure data were available. A 
cross-plot of calorific value (MMMF) versus the calculated 
depth of burial (Fig. 13) follows Hilt's rule, namely, that 
with an increase in depth of burial, coal increases in rank. 
Hacquebard (1977) indicated a logarithmic relationship 
between equilibrium moisture and original depth of burial 
(Fig. 12, Equ. 3), and the present study documents an 
inverse linear relationship between equilibrium moisture and 
calorific value (MMMF). Therefore, the relationship between 
Original depth of burial and calorific value should not be 
linear. However, because of the small range of depth of 
buriele values calenlatedr (4 200°to 2200°m (3937: to 7218 ft]), 
the relationship appears linear, (especially for the 
calorific range 22000 to 29000 kJ/kg) and a simple least 
Squares regression analysis can be performed. Results of 


this regression analysis yields the equation 
CV(MMMF) = 14748 + 6.25(DOB) (4). 


The slope component of this equation indicates that plains 
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Figure 13: Crossplot of calorific value (moist mineral 
matter free basis) versus original depth of burial. 
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coals increase in calorific value (MMMF) at a rate of 

6.25 kJ/kg per metre (0.82 Btu/lb/ft) of depth (or 
overburden). This rate of increase represents the 
coalification gradient. Equation (4) or the coalification 
gradient should not be used for coals with ranks beyond the 
calorific value range of 22000 to 29000 kJ/kg because of the 
curvilinear relationship between depth of burial and coal 
rank. 

Campbell (1972), using coal samples from deeper oil 
wells, calculated a coalification gradient of 5.34 kJ/kg per 
metre (0.7 Btu/lb/ft) for the Alberta plains area. This 
lower gradient, however, was calculated on the basis of an 
assumed linear relationship between rank and depth, for the 
range of coal ranks between subbituminous and low volatile 
bituminous. As discussed previously, however, the assumption 
of a linear relationship for such a considerable calorific 
value range cannot be made. 

With the use of Equation (4), the overburden required 
to produce various low rank coals can be calculated 
(Table 11). This equation can also be used to modify 
Hacquebard's (1977) diagram (Fig. 12) to include the coal 
ranks found in the Alberta plains (Fig. 14). 

It should be noted that in order for this relationship 
to hold true, sediment deposition between the oldest and 
youngest coals discussed (Belly River Group and Scollard 
Member, respectively) would have to have been generally 


continuous, without major periods of erosion. The effect of 
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Figure 14: Modification of Hacquebard's (1977) diagram, 
including the coal ranks found in the plains area of 
Alberta. 
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TABLE 11 - Required thickness of overburden to mature coals. 
The calculation for subbituminous A coal assumes that the 
coal is non-agglomerating. 


Coal Rank (ASTM) Required Overburden 
(metres) 

Subbituminous C =o crtesgs Sea A 

Subbituminous B 11760 =) 1548 

Subbituminous A $5426e— 1920 

HV Bituminous C TO 20a ums > 


a major period of deposition, and subsequent erosion, would 
be that the coals underlying the erosion surface would have 
anomalously higher ranks, as compared to coals overlying the 
erosion surface. This would occur because the coals from 
below and above the erosion surface would have had different 
compaction histories, and thus reached their peak maturities 
at different times. Major unconformities, which would 
reflect periods of significant erosion, have not been 
recognized within the stratigraphic section being discussed. 
Therefore, in a vertical section (for example, a borehole), 
it is assumed that all coals reached their maximum (and 
current) level of maturation at the same time; this time 
being at a point when their depth of burial was greatest. 
Once the original depth of burial of a coal sample has 
been calculated, the amount of sediment removed above the 
coal sample can be determined by subtracting the present 


depth of the sample from its original depth of burial. In 
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order to present this data in map view, and to remove some 
of the statistical variation in the data, the coal-bearing 
portion of the plains (as defined by Fig. 1) was divided 
into blocks, each nine townships in area (840 km? [324 
mi?]), and an average value of removed overburden calculated 
for each block (Appendix 4). These values were then plotted 
and contoured (Fig. 15). 

As is shown in Figure 15, the amount of overburden 
removed varies from approximately 900 m (2,953 ft) in the 
east to approximately 1900 m (6,234 ft) in the west part of 
the plains area. 

Suggate and Lowery (1982) note that coals with similar 
equilibrium moisture contents have higher ranks in areas 
with higher geothermal gradients compared to coals from 
areas with lower geothermal gradients. As discussed in this 
thesis, coals with similar equilibrium moistures must have 
been subjected to similar compaction histories (or similar 
depths of burial). The reason, therefore, that coals from 
areas with higher geothermal gradients are higher in rank 
would be because these coals reached higher levels of 
chemical maturation. Comparing figures 12 and 14, the coals 
occurring in the plains area of Alberta (Fig. 14) appear to 
mature at shallower depths, as compared to the coals used in 
the Hacquebard (1977) study (Fig. 12). This Pood indicate 
that the regional geothermal gradient, at the time of 


maximum burial, was possibly higher in the plains area of 
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Alberta. ' 

The high volatile bituminous C coals in this study were 
subjected to an estimated paleotemperature of between 107 
and 99 degrees C, using models by Bostick (1979) and Price 
(1983). Assuming a surface temperature of 16 degrees C for 
the time the maximum paleotopographic surface existed 
(Hacquebard, 1977), and using a depth of burial of 1920 m 
(6299 ft) for coals that have an estimated vitrinite 
reflectance of 0.51%, these temperatures translate to a 
paleogeothermal gradient of between 47.3 and 43.2 degrees C 
per km of depth. These gradients are considerably higher 
than would be expected. Possible reasons for these high 
estimates are: 1) the models provided by Bostick (1979) and 
Price (1983) are incorrect, or do not apply to the plains 
area of Alberta, 2) the estimate of vitrinite reflectance 
from Berkowitz (1979) is too high for the coals discussed, 
or, 3) the thicknesses of overburden estimated are eee than 
the actual thicknesses that existed. Further work on the 
reflectance of plains coals is required in order to provide 
more detailed information on paleogeothermal gradients. 

Cross-section AA' (Fig. 16) is constructed through the 
plains, from northeast to southwest (Fig. 1), and indicates 


the subsurface stratigraphy, the westerly dip of the 


‘For reference, the present geothermal gradient in the 
plains area of Alberta decreases from east to west, from 
approximately 33 degrees C to 22 degrees C per km (United 
States Geological Survey, 1976), and Hitchon (in press) 
Suggests that this geothermal gradient pattern existed since 
early Eocene time. 
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coal-bearing units, and Pheu peconscructed maximum 
paleotopographic surface. The maximum paleotopographic 
surface was determined by adding the calculated amount of 
removed overburden to the present land surface. During the 
period of time when the calculated maximum paleotopographic 
surface existed, plains coals would have reached their 
maximum (and current) degree of metamorphism. As noted in 
Table 11, high volatile bituminous C coals require at least 
1920 m (6299 ft) of burial. A line 1920 m (6299 ft) below 
the maximum paleotopographic surface, shown in cross-section 
AA' (Fig. 16), should therefore separate the occurrence of 
bituminous coals from subbituminous coals in the subsurface. 
A similar method of calculation can be used to determine the 
theoretical occurrence of the subbituminous ranks of coal in 
the subsurface. This has been done in cross-section BB', 
which is constructed east to west through southern Alberta 
(Fig. 17). As noted in both cross-sections AA' and BB' 
(Figs. 16 and 17), higher rank coals occur nearer the 
present land surface in the western part of the plains where 
calculations show that greater amounts of sediment have been 
removed. 

At first glance, the thicknesses of removed overburden, 
as noted in Figures 15, 16 and 17, appear excessive. 
However, the following discussion on Tertiary he aa 
and tectonic events may demonstrate that these thicknesses 


are realistic. 
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Taylor et al]. (1964) indicate that Paleocene sediments 
(lower Tertiary) attain a thickness of over 1500 m (4921 ft) 
in portions of the westernmost interior plains and it is 
very possible that equivalent sediments extended eastward in 
the past. In addition to these sediments, sediments derived 
from the uplift of the Rocky Mountains during the Eocene and 
Oligocene Epochs (Taylor et a/]., 1964), would have provided 
additional sediment thickness. Taylor et al]. (1964, p. 193) 
state 


",..-only an insignificant portion of the enormous 
amount of debris that moved eastward along drainage 
channels in Eocene and Oligocene times remains on 
the Canadian Plains". 

Price et al. (1970) differ from Taylor et a]. (1964) in 
their interpretation of the structural development of the 
eastern Rocky Mountains. They feel that the major thrusting 
culminated during the Paleocene, *not*during the Eocene and 
Oligocene as suggested by Taylor et al. (1964). The 
deposition of Eocene - Oligocene sediments in the plains 
area was, by Price et a]. (1970, p. 23) 

",.,.associated with a new and different phase of 

deformation in the Cordillera - gravity faulting, 

during which block “fault *structures were «superposed 

over older thrust and fold structures." | 

Regardless of their views on the structural development 
of the easternmost Rocky Mountains and westernmost plains, 


both authors support the idea of a blanket of Eocene - 
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Oligocene sediments covering the plains area. Significant 
portions of this Tertiary sediment, as well as portions of 
Upper Gretaceous sediment, were subsequently eroded, as 
shown by the present elevation and stratigraphy of the 
Cypress Hills (Fig. 1). 

The east portion of cross-section BB' (Fig. 17) 
displays both the present surface elevation through township 
13 and the elevation of the Cypress Hills, through 
township 8 (cross-section CC', 50 km [31 mi] south of 
cross-section BB'). The maximum paleotopographic surface in 
cross-section BB', projected easterly, would fall 
approximately 180 m (591 ft) above the top of the Cypress 
Hills. The Cypress Hills are capped with a coarse Oligocene 
conglomerate, and it is possible that at least 180 m 
(591 £t) of sediment could have been eroded from the top of 
the these hills during the last 30 million years. 

Since the calculations of sediment removal in the 
vicinity of the Cypress Hills appear reasonable, the 
estimated sediment removed in the contoured area is probably 
reasonable (Fig. 15). 

Geophysical logs from petroleum boreholes located in 
the Cypress Hills (6-3-8-2-W4M), and outside the hills 
(11-16-13-3-W4M) on cross-section BB', were examined to 
determine the thickness of rock eroded in areas re from 
the Cypress Hills (Fig. 18). Comparing the two boreholes, 
and using the Geological Map of Alberta, indicates that 45 m 


(148 £t) of Belly River Group, 336 m (1104 ft) of Bearpaw 
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Figure 18: Cross-section DD' (see text for cross-section 
location), indicating the thickness of units eroded 


in the vicinity of the Cypress Hills. 
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Formation, 57 m (186 ft) Oieueasechce Canyon Formation 
equivalent', 7 m (22 ft) of Battle Formation, and at least 
218 m (717 £t) of post-Battle Formation sediments have been 
removed from the vicinity of the 11-16 borehole. Because the 
above mentioned lithologic units are only exposed ina few 
wells in the Cypress Hills area, trends of formation 
thickening or thinning could not be determined, and 
calculations of sediment removal assume that all rock units 
remain constant in thickness. Excluding the thickness of 
sediment that has been removed since Oligocene time from the 
Cypress Hills, the calculations show that at least 678 m of 
sediment have been removed from the area of the 11-16 
borehole. 

Calculation of meaningful denudation rates require 
knowledge of the area of the source area, the volume of 
sediments derived from the source area, and the duration of 
the erosion (Menard, 1961). The thickness of eroded sediment 
calculations in this study cover only a portion of a 
drainage basin, thus the volume of sediments removed from 
the drainage basin cannot be determined. Menard (1961) 
calculated a denudation rate of 46 mm (1.8 in.) per 1000 
years for post-Jurassic erosion of the Mississippi drainage 
basin. Accurate paleodraingage directions for Oligocene and 
Eocene sediments in the plains area are not nave because of 
the scarcity of outcrop, however, early Tertiary drainage 
did have a east to southeast component (Rahmani and 


‘known as the Ravenscrag Formation in the Cypress Hills 
area, following Carrigy (1971) 
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Lerbekmo, 1975). If these Haauenal paleodrainage conditions 
persisted through middle Tertiary time, a portion of the 
sediments indicated as eroded in this study would have 
entered the Mississippi drainage basin. Assuming an erosion 
time of 30 million years (middle Oligocene to present), the 
central plains area of Alberta underwent erosion at rates 
between 33 and 63 mm (1.3 and 2.5 in.) per 1000 years. These 
values match well with the overall denudation rates 
calculated for the Mississippi drainage basin. 

Magara (1978), using sonic transit times in shales to 
estimate the degree of compaction, also calculated the 
thickness of sediment that was eroded from the plains of 
Alberta. The thickness he calculated, however, is 
considerably less than the thicknesses presented here 
(Figs. 15, 16 and 17). Magara's (1978, p. 44) cross-section 
III-III' trends approximately northeast in the southern part 
of Alberta. At a point on this cross-section, approximately 
14 km (9 mi) north of the city of Medicine Hat (point M in 
Fig. 1), Magara (1978) calculated approximately 176 m 
(577 £t) of sediment erosion. The present surface elevation 
atithis point hish732 m (2,400 8t) hhindicatiing ‘a maximum 
paleotopographic elevation of 908 m (2,977 ft). The Cypress 
Hills, 68 km (42 mi) southeast of point M (Fig. '), and 
approximately on depositional strike, have an elevation of 
1465 m (4,807 ft). These two locations should have had 
Similar elevations during Oligocene time when sediments were 


actively being shed easterly from the foothills and 
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mountains. Magara's (1978) estimation of maximum 
paleotopographic elevation at point M may be at least 557 m 
(1,827 ft) less than the possible maximum, as suggested by 
the present day elevation of the Cypress Hills. 

Because of this discrepancy, the remainder of Magara's 
calculations of sediment erosion in the Alberta plains may 
be equally at variance with those of this study. 

Coals equivalent in age to those occurring in the 
Scollard Member in the plains area are being mined in the 
Coal Valley area in the foothills (Fig. 1). Equilibrium 
moisture content of these coals iS approximately 8.5 percent 
and calorific values (MMMF) vary from 27905 to 28424 kJ/kg 
(11,997 to 12,220 Btu/lb) (Mervyn Rogan, Luscar Ltd., per. 
comm., 1982). Coal quality parameters approaching these 
values can be noted in coal samples from Upper Cretaceous 
and Tertiary coals in the plains (Tables 2 and 7). If, as 
previously thought (Campbell, 1983; Stansfield and Lang, 
1944), coal rank in the plains area of Alberta was related 
to heat and / or pressure generated from the uplift of the 
Rocky Mountains, coals in the foothills region should be 
considerably higher in rank than those found in the plains. 
This, as shown above, is not the case. 

The similarity between coals found in the Coal Valley 
area and those found in the subsurface of the central plains 
area can be explained by the model for coal rank 
distribution presented in this study. The Coal Valley area 


was uplifted during Paleocene or early Eocene time (Price et 
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al., 1970; Taylor et a/].,1964) and thus coals in this area 
would not have had the maturation benefit of the overburden 
added during Eocene and Oligocene time. Upper Cretaceous and 
Tertiary coals, therefore, increase in rank in the plains 
from east to west because of increasing original overburden 
thicknesses, and then decrease in rank when encountered in 
the foothills region, because of thinner, or non-existent, 
post-Paleocene age overburden thicknesses. 

In a previous section, where the relationship between 
equilibrium moisture and ash was discussed, it was noted 
that in a study of Saskatchewan lignites, equilibrium 
moisture content and ash content had a correlation 
coefficient of R = -0.90, considerably higher than the 
correlation found in this study (where R = -0.38). This 
study demonstrated that the variation in equilibrium 
moisture in coals from the plains area of Alberta can be 
related to the variation in coal seam depth of burial,-: which 
varies according to the sample location. Because of the good 
correlation between equilibrium moisture and ash in the 
Saskatchewan study, other factors affecting the lignites 
equilibrium moisture, such as variation in depth of burial, 


did not exist, or played a lesser role. 
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V. ECONOMIC IMPLICATIONS 

The relationship between calorific value (MMMF) and 
depth of burial (Equation 4) was used to determine the 
thickness of overburden required to mature various low rank 
coals (Table 12). This information, combined with the map 
showing the thickness of removed sediment (Fig. 15), can be 
used both to define areas containing similar ranks of coal 
and to determine present depths to coals of higher ranks. As 
high volatile bituminous C coals require at least 1920 m 
(6,299 £t) of burial (Table 12), these coals should be found 
at or near the surface in areas where approximately 1920 m 
(6,299 ft) of overburden have been removed. Similarly, these 
coals should be found approximately 100 m (328 ft) below the 
surface where 1820 m (5,971 ft) of overburden have been 
removed. Using this method we can calculate the minimum 
theoretical depth to the occurrence of high volatile 
, bituminous C coal in the plains (Fig. 19) where near-surface 
data exists. Although the maps are not presented here, a 
Similar method of calculation can be used to define the 
occurrence of the subbituminous ranks of coal. 

According to ERCB (1982), 82% of the total area 
underlain by coal-bearing strata in Alberta occurs in the 
plains area south of Lesser Slave Lake. Indications from 
maps and cross-sections (Figs. 15, 16 and 17) in this study 
are that the thickness of eroded sediment increases in a 
westerly direction, therefore coals occurring in that 


direction would increase in rank. This is in part shown by 
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Figure 19: Calculated depth to high volatile bituminous 
C coals in the plains area of Alberta. 
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the higher average rank of Scollard Member coals versus 
Horseshoe Canyon Formation coals (Table 7; subbituminous A 
versus subbituminous B, respectively), with Horseshoe Canyon 
Formation coals generally sampled to the east of Scollard 


Member coals. 
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VI. CONCLUSIONS 

Upper Cretaceous and Tertiary near-surface coals in the 
Alberta plains vary in rank from subbituminous C to high 
volatile bituminous C. Considering coals that contain less 
than 25% ash (D), equilibrium moisture averages 17.8%, ash 
(D) averages 14.4%, volatile matter (DAF ) averages 43.0%, 
fixed carbon (DAF) averages 57.0%, oe Pee value 
(MMMF) averages 23967 kJ/kg. Calorific value is highly 
dependent on ash and moisture content. 

Coal rank in the Alberta plains increases with depth, 
at a rate of 6.25 kJ/kg per m (0.82 BTU/1b/ft) and increases 
towards the Rocky Mountains. The former increase is related 
to the coalification gradient and the latter to greater 
depths of burial. Subsequent to the coals reaching their 
maximum maturity in the mid Tertiary, erosion has removed 
between 900 and 1900 m (2953 and 6234 ft) of sediment from 
the plains area. Erosional outliers, such as the Cypress 
Hills, appear to fit the proposed erosion model. 

The results of this study suggest that a significant 
portion of Alberta's Upper Cretaceous and Tertiary coal 
resources, which traditionally have been classed as 
subbituminous in rank, may fall into the bituminous ranks. 

The hypothesis proposed in this study could be tested 
by sampling coals from a series of coal testholes or o11 and 
gas wells in an east - west direction. A combination of 
equilibrium moisture, vitrinite reflectance and rank data 


from these samples would provide information on 1) the 
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thickness of removed strata both east and west of the 
present study area, 2) information on paleogeothermal 
gradients, and 3) the rate of rank increase with respect to 


both depth and geographic location. 
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A. APPENDIX 1 


COAL ANALYSIS FILE 


The 685 analysis enclosed in Appendix 1 are sorted 
firstly by land location, from south to north, and secondly 
by depth. The land locations on the following pages are 


presented in their conventional format, namely 


LSD-SEC-TWP-RGE-MER 


In order to minimize the amount of space required for 


Appendix 1, the following abbreviations are used. 


Abbreviation units Parameter 
COAL ZONE Coal-bearing Unit 
SC-M Scollard Member 
HC-F Horseshoe Canyon Formation 
BR-G Belly River Group 
WP-F Wapiti Formation 
UNKN Unknown unit 
UPPER DEPTH m upper sampling depth 
LOWER DEPTH m lower sampling depth 


SAMP TYPE sample type 
CUTT cuttings 
CORE core 
MOISTURE moisture 
AD % as determined moisture 


AR 
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as received moisture 
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MEQ % equilibrium moisture 
ASH DRY % dry basis 
VL-MATTER volitile matter 
DRY % dry basis 
DMMF % dry mineral matter free basis 
F X-CARBON fixed carbon 
DRY % dry basis 
DMMF % dry mineral matter free basis 
CAL-VALUE calorific value 
DRY kJ/kg dry basis 
MMMF kJ/kg moist mineral matter free basis 
ULT-ANALYSIS ultimate analysis 
S % sulfur, dry basis 
C % carbon, dry basis 
H % hydrogen, dry basis 
N % nitrogen, dry basis 
O % oxygen, dry basis 
COAL RANK coal rank determinations 
HVB-C high volatile bituminous C 
SUB-A Subbituminous A 
SUB-B Subbituminous B 
SUB-C Subbituminous C 


INSF insufficient data available 
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B. APPENDIX 2 


given base 


as received 
as received 


as received 


dry 
dry 
dry 
Where: 
MAR 
MEQ 
AAR 
AD 


is 


CONVERSION CHART 


required base conversion 
dry 100/(100-MAR) 
dry ash free 100/(100-MAR-AAR) 
moist (100-MEQ) /(100-MAR) 
as received (100-MAR) /100 
dry ash free 100/(100-AD) 
moist (100-MEQ) /100 


= as received moisture 
= equilibrium moisture 


= ash, as received basis 


= ash, dry basis 


Multiply the analysis component to be converted by the 


conversion factor in order to obtain the required base. 


ATHAM-OOT) OGL 
O00 f )% (OSM-G0T) 


Ot \ (HfAMHDO Ph) 


Gam-001 )\00T, © 


f 


7. 
7 
7 


J 


at 


nelerevaos 


») 
i 


a 


na 


GA-OOF )\ 90 


t\ (OSM 00T ) 


sivjelom bavigesy 66 
sivdeiom miltmit rags 
elesa bev isoss es «tae 


2iesd ¥rbh ea 


yd Ssirsvnos-sd 03 shilpa sisq inne sad i yiqivtun 


2286: Peainpsx eds obsiti ov rabao ai tesa} solevevnos | 


ae 


116 


C. APPENDIX 3 


PARR FORMULAE 


DMMFFC = [ (Fc-0.15S)(100)]/[100-(MEQ+1.08A+0.55S) ] 

DMMFVM = 100-DMMFFC 

MMMFCV = [ (cv-50S) (100)]/[100-(1.08A+0.55S) ] 
Where: 


DMMFFC = dry mineral matter free fixed carbon 
DMMFVM = dry mineral matter free volatile matter 
MMMFCV = moist mineral matter free calorific value, BTU/1b 
FC = fixed carbon, moist basis 
S = sulfur, moist basis 
MEQ = equilibrium moisture 
A = ash, moist basis 


CV = calorific value, moist basis, BTU/1lb 


To convert BTU/l1b values to kJ/kg, multiply the BTU/1b 


value by a conversion factor of 2.326. 
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D. APPENDIX 4 


REMOVED OVERBURDEN DATA 


The location given (meridian, township, range) 
represents the centre of a 9 township block. The removed 
overburden value represents the average removed overburden 


for the 9 township block. 


Mer Twp Rge Removed overburden (m) 
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